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O . Abstract 

■ We construct effective 3d field theories for the Minimal Supersymmetric Standard 

p • Model, relevant for the thermodynamics of the cosmological electroweak phase transi- 

p . tion. The effective theories include a 3d theory for the bosonic sector of the original 

D ! 4d theory; a 3d two Higgs doublet model; and a 3d SU(2)+Higgs model. The inte- 

^ ' grations are made at 1-loop level. In integrals related to vacuum renormalization we 

take into account only quarks and squarks of the third generation. Using existing non- 
^ ■ perturbative lattice results for the 3d SU(2)+Higgs model, we then derive infrared safe 

upper bounds for the lightest Higgs boson mass required for successful baryogenesis at 
the electroweak scale. The Higgs mass bounds turn out to be close to those previously 
found with the effective potential, allowing baryogenesis if the right-handed stop mass 
parameter mfj is small. Finally we discuss the effective theory relevant for very 
small, the most favourable case for baryogenesis. 
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1 Introduction 



The generation of the baryon number of the Universe remains to be satisfactorily 
explained. It is quite plausible, though, that an important role in the process was played 
by the cosmological electroweak phase transition Within the Standard Model the 
phase transition appears nevertheless to be too weakly of first order to produce the 
baryon asymmetry for realistic Higgs masses (for a review, see 0). Additional problems 
may be related to the amount of CP-violation available. Hence one is led to extensions 
of the Standard Model. 

One possible consistent extension of the Standard Model is the Minimal Supersym- 
metric Standard Model (MSSM). It has a large parameter space available so that it 
should be possible to find some corner with a strong enough first order transition. In 
addition, there are additional sources of CP-violation. Indeed, the electroweak phase 
transition in MSSM has been studied quite actively [3-7]0. 

The investigations made so far (apart from I^H], ^) have been based on the 1- and 
2-loop effective potentials for the Higgs field. The limit that the CP-odd Higgs mass 
rriji is infinite was taken in [Q, ^, ^ |^, leaving just one Higgs doublet and being the 
most favourable case for baryogenesis |Q. The result of these investigations was that 
in general, it appears difficult to make a strong enough transition unless the right- 
handed soft supersymmetry breaking stop mass parameter mfj is small. Recently it 
has been noted that even smaller values of mfj than originally considered should be 
phenomenologically possible , leading to a transition that is definitely strong enough 
for baryogenesis. In addition, 2-loop effects have been found to be favourable 0. 

All the studies based on the effective potential are subject to the infrared (IR) 
problem at finite temperature [§]. The IR problem is related to the zero Matsubara 
components of bosonic fields, and precisely these components account for the cubic 

1- loop terms in the effective potential studied in [3-6], as well as for the logarithmic 

2- loop terms making the effect in [0 . The IR-problem calls for non-perturbative inves- 
tigations of the problem. The method of choice for non-perturbative investigations is 
the framework of dimensional reduction [9-17]. Dimensional reduction means that one 
constructs an effective 3d theory producing the same Green's functions as the original 
theory for the light bosonic fields. The perturbative dimensional reduction step is free 
of IR-problems, and the resulting super-renormalizable 3d theory can then be studied 
with high precision Monte Carlo simulations [18-22]. 

The non-perturbative investigations of the electroweak phase transition in the Stan- 



dard Model have revealed the following pattern [IS]. As long as the transition is strong 



enough for baryogenesis, the IR-problems are not very dramatic and effective poten- 
tial studies do produce a reasonable estimate of the properties of the phase transition. 



^ Upon completion of this work, three more papers on the same subject appeared [37-39]. In 
the 1-loop effective potential is studied. In ||3|, ^ the authors study the dimensional reduction of 
MSSM as in the present work. In |Q the analysis is a bit less complete than here and the conclusions 
are somewhat different. In ]39[ | the formulas for dimensional reduction and heavy scale integrations 
are in some parts more, in some parts less complete than here, but vacuum renormalization and the 
implications of the formulas to the electroweak phase transition are not discussed. 
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When the transition gets weaker, non-perturbative effects become large. However, even 
if non-perturbative effects are small for stronger transitions, it is nevertheless interest- 



ing to note that prior to the lattice study in [|I9| many perturbative studies stated 
that baryogenesis was possible up to mn ~ 45 GeV. In it was discovered that 
practically no Higgs mass is possible. While this effect is mostly related to vacuum 
renormalization instead of non-perturbative IR-effects, it nevertheless proves that it is 
important to work in a consistent framework where all the approximations are under 
control. 

The purpose of the present paper is to make a dimensional reduction for the MSSM. 
We also perform further integrations inside the dimensionally reduced 3d theory, to 
arrive at the simplest possible effective theory. In particular, we construct a 3d two 
Higgs doublet model and a 3d SU(2)+Higgs model in the part of the parameter space 
where it is possible. For the latter theory, the existing non-perturbative lattice results 
allow to remove the IR problem from the Higgs mass bound. The bound derived is 
in principle also gauge and /2- independent, unlike the ratio v{Tc)/Tc derived from the 
effective potential. 

On the technical side, one purpose of the present investigation is to study how the 
cubic scalar vertices, not present in the Standard Model, affect dimensional reduction. 

In comparison with [4-7] , we also try to be more explicit about the effects of vacuum 
renormalization. The theory studied is more or less the same. We include here the 
bottom Yukawa coupling hb and study a general CP-odd Higgs mass niA as in 

It is found that in the region of the parameter space where reduction into the 3d 
SU(2)-|-Higgs model is possible and the transition is strong enough for baryogenesis, 
the non-perturbative results agree with the effective potential investigations. The con- 
clusion is that Higgs masses mh<75 GeV produce a strong enough transition if mfj 
is small enough, m^<(50 — 100)^ GeV^. Hence, the situation has improved with re- 
spect to the Standard Model where no Higgs mass is possible. Where reduction into 
SU(2)+Higgs cannot be made — notably when is still smaller and the transition is 
even stronger — we propose an effective 3d theory allowing more detailed studies 
of the problem. 

The plan of the paper is the following. In Sec. ^ we briefly review the Higgs mass 
bound in the Standard Model and its derivation within the 3d framework. In Sec. |^ 
we state in some detail the approximations adopted and the Lagrangian used in the 
present investigation. Sec. ^ contains the dimensional reduction into a 3d bosonic 
effective theory. In Sec. ^ we make further integrations inside the 3d theory, removing 
the squarks and the temporal components of the gauge fields. The resulting two Higgs 
doublet model is diagonalized in Sec. |], and the heavy Higgs doublet in integrated out 
in Sec. 0. In Sec. ^we discuss how the running Lagrangian parameters are fixed through 
vacuum renormalization. The numerical results for the strength of the transition are 
in Sec. ^. Finally, in Sec. |TD| we propose an effective theory for describing the phase 
transition if the mass parameter mfj is very small. Sec. |ll] is the conclusions. 
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2 The EW phase transition in the Standard Model 



The thermodynamics of the electroweak phase transition in the full Standard Model 
has been extensively studied in the literature (0 and references therein). Perturbative 



studies exist up to 2- loop level p3|, ^ Non-perturbative lattice studies rely 



on perturbative 2-loop dimensional reduction [12-15], and have been performed for a 
wide range of Higgs masses [18-22]. The Higgs mass bound in terms of the parameters 
of the 3d SU(2)+Higgs model was derived in [IS . 

The Higgs mass bound arises as follows. Assume that there in some underlying 
physical 4d theory in which the electroweak phase transition takes place so that the 
static Green's functions of the lightest excitations are described by the effective theory^ 

L3d = ^i^^.i^^. + (A0)^(A0) + m20V + A3(0V)', (2.1) 

where Di = di—ig3T'^A^/2. Then the phase transition is strong enough for baryogenesis 
if at the phase transition point [19 



X = ^ < 0.03 - 0.04. (2.2) 

93 

Since the theory in eq. (|2.1|) is super-renormalizable, the parameters A3, (^f do not run 
and the quantity x is a well-defined pure number. It is also gauge- independent. The 
uncertainty in ( p.2|) arises from uncertainties in estimates of the sphaleron rate in the 
broken phase and from uncertainties in the real-time dynamics of the phase transition 
(whether the Universe reheats back to Tc after the nucleation period, etc.). 
For the Standard Model, the parameters m|, A3 and g"^ have been calculated in terms 



of temperature and the physical zero-temperature parameters of the theory in [[15 
Then one may solve for the critical temperature from the condition 

ml{mh,T,) = 0, (2.3) 

and use this in the estimate of 

_ A3(m,,r,) 

^5i(m.,T,)- ^'-^^ 

Eq. ( |2.3| ) does not give Tc exactly (it corresponds to resummed 1-loop accuracy), but 
this does not matter since Xs/g^ depends on Tc only through logarithmic 1-loop cor- 
rections. From an analysis of the type outlined, one gets that no Higgs mass (or at 
most an extremely light Higgs mass, m/j<20 GeV) would satisfy the bound ( |2.2| ) in the 
Standard Model since x > 0.04 due to top Yukawa coupling corrections, see Fig. 27 

m. 



m 



In this paper we study whether a theory of the type in eq. (|2.1| ) can be constructed 
in the MSSM and what would be the Higgs mass bound implied. 



^ It should be noted that even though only the SU(2) group is displayed explicitly in eq. {2A), the 
perturbative effects of the U(l) group, making the phase transition stronger, have been included in 
the bound ( |2.2| ). No non-perturbative lattice simulations exist yet for the SU(2)xU(l)-|-Higgs theory. 
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3 The Lagrangian 



We start by discussing the Lagrangian used and the simphfications made, fixing at the 
same time the notation. We work throughout in Euchdian space and for definiteness 
in the Landau gauge. The value of x derived is gauge- independent. 

The main simphfications are the following (for the complete Lagrangian in MSSM, 
see [^). First, we neglect the U(l) subgroup in loop corrections related to vacuum 
renormalization and dimensional reduction. That is, no difference is made between 
and + g'"^ beyond tree-level in IR-safe integrals. This is a good approximation 
as far as the electroweak phase transition is concerned, especially with respect to the 
other uncertainties in the calculation. Even at tree-level, we display explicitly only the 
covariant derivatives related to SU(2) and SU(3). 

Second, we will assume that the gaugino and higgsino mass parameters in the sym- 
metric phase are so large that these fields have decoupled, as is usually assumed in the 
present context [4-7]. Even if the masses are smaller, these fields do not have very 
much significance, being fermions: at finite temperature, the important effects arise 
from IR-sensitive bosons. In the framework of the present paper, the extra fermions 
would only affect the parameters of the 3d theory in the first dimensional reduction 
step, but the later integrations remain precisely the same. It should also be noted that 
gauginos and higgsinos do not couple to the scalar Higgs degrees of freedom through 
the dominant Yukawa coupling /i^, unlike the top quark. In general, it is expected that 
the effect of gauginos and higgsinos would be to make the phase transition weaker, 
due to the increased screening in the thermal masses ^, 0]. However, gauginos and 
higgsinos do have an effect when 1-loop corrections to the top Yukawa coupling ht are 
calculated. Since ht gives the most important effects in the present calculation, the 
loop corrections may also be important. We return to this point in more detail below. 

Third, only the squark partners of top and bottom quarks are assumed to be light 
enough to affect the electroweak phase transition. 

Then the remaining fields are as follows: there are the SU(2) and SU(3) gauge fields 
A^, C^. The Higgs fields are H^,H^, with hypercharges Y = — The adjoint 
Higgs fields with opposite hypercharges are denoted by 

= iT2H^\ H'^ = iT2H'^*. (3.1) 

The 2-index antisymmetric tensor is defined through ei2 = —1, so that = —tijH^*. 
We use the notation 

ff^-(Hl\_ 1 / hl + ^hl\ rr2_(Hf\_ 1 (hl+thi\ 

~[m )~ V2[ -hl + ^hl j ' )- ^[ hl-^hl ) ^^"'^ 

for the complex and real components of the Higgs fields, so that at zero temperature 

{hi) = v„ (hi) = V2. (3.3) 
The fermions of the third generation are 



QLa — T^" , tRa, bRa, (3.4) 
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where a is the SU(3)-index and the hypercharges are Y = 1/3,4/3, —2/3, respectively. 
Correspondingly, the squarks of the third generation are 



Qc 



ha 



t 



(3.5) 



with the hypercharges Y = 1/3,-4/3,2/3. The fields U,D transform under SU(3) 
with the adjoint generators = — A^. 

The part of the action containing the kinetic terms of and interactions between gauge 
fields and fermions remains as in the Standard Model. For the Higgs and squark fields 
the quadratic terms are 

L = iD;HyiD;H') + iD;HyiD-H') 

+ (D-Q)t(2;-Q) + {DIW)\D^JJ*) + {DID*)\D^^D*) 



(3.6) 



where D"^,^ 



dfj, — igr"-A1j^/2 — and w and s indicate the charge included. 

The supersymmetric interactions are generated by the superpotential and by the 
Z)-terms. We take the superpotential to be 



W = ^eijHlHj + hteijHfQjaUa + hbeijU^QjaDa- 



(3.7) 



Hence also the bottom Yukawa coupling hf, is kept, although its effect is small since the 
region of parameter space which can affect baryogenesis is around tan/3 = V2/V1 ~ 2, 
and hb ~ nib / {2m z cos f3). The interaction Lagrangian following from the superpoten- 
tial is 



Lw = ht (tRH^^qL + qiEhR) +h(bRH^^qL + qLH^hR^ 
The interaction Lagrangian following from the D-terms, on the other hand, is 



(3.8) 



+ 
+ 



/O* /O I Tjm* Tjm I Tjn* rrn 



'ij'kl 



9l 



Here we kept the U(l) coupling g' only in the Higgs sector. In eq. 

A A 2 

'Tij'Tkl = C^^il^jk - Sij6kl), X^ijX^s = -{SSasSp^ - Sa/sS-ys)- 



(3.9) 



(3.10) 



The soft supersymmetry breaking cubic interactions are 



u 



>H^^QaUa + d.H'^Q^D^ + CsH^^QcD^ + w^H'^Qo^Uo^ + H.c. (3.1i; 



Here we use the notation of with opposite signs. The relation to the more standard 
A-parameters is discussed in Sec. || in connection with vacuum renormahzation. In 
particular, it should be noted that the parameter includes the term fihb arising 
from the superpotential and the parameter Wg includes —^ht. Since there is also an 



arbitrary soft component in these terms |26| and since we assume the higgsinos to be so 



heavy that they have decoupled, the theory does in fact not depend at all on the true 
super symmetric mass parameter /i. In the following, we shall restrict the parameters 
Us,ds,es, Ws to be real. 

Although quite a few simplifications have been made, there are still a lot more 
parameters left than in the Standard Model. The two scalar sector parameters i>^,X 
appearing there are replaced by , m^, ^^2, rriq, mfj, rrijj, Us,ds,es,Ws- 

There is the following important point to be noticed about the coupling constants 
g, gs, ht, hb in the present theory. As an example, take the gauge coupling. If one takes 
the theory under investigation as such, then the weak gauge coupling in the gauge 
sector runs as 

_ dg\fi) g'' SriF + N^-AA 

^ dfi ~8n^ 6 ' ^^-^^^ 

where n^r = 3 is the number of families and A^<j = 5 is the number of scalar doublets 
interacting with the SU(2) gauge fields. However, the couplings gf, g2, di, dl = g'^ in 
the SU(2)-part of the scalar potential following from eq. (|3.9| ), 



V = -gl{H^^H^f + -gl{H^^H^f + -g^H'^ H^H^^ - - gjH^^ H'^ H\ (3.13) 

8 8 4 2 

run as 

- dgfifi) ^ - dg^ifi) _ _ dgl{fi) _ 7 g^ 

dji dfl dfi 2 87r^ ' 

- dglifi) 5 g^ 

= -28^- ^^-''^ 

Hence within the present theory one would have to renormalize these couplings sep- 
arately from the coupling in the gauge sector. In other words, one has to consider a 
large number of zero-temperature observables in terms of which to fix the independent 
parameters. If on the other hand one wants to maintain the universality of the gauge 
coupling, then one has to include the complete supersymmetric structure of the theory 
in the calculation in one way or the other. In the present theory, supersymmetry is 
maintained only in the quark-squark sector of the third generation, and indeed, if only 
these fields are included in the internal lines of loop integrals, then g"^ runs everywhere 
as 

-dg^fi) Sg' 

= 2 8^ ^'-''^ 
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We shall work within the accuracy of this approximation here and assume the gauge 
coupling to be universal. 

The same thing applies also to the Yukawa couplings ht, hf, and is quite important 
there as well, since ht is large and gives the dominant effects. Indeed, within the present 
theory, the Yukawa couplings in the different squark-Higgs and quark-Higgs interactions 
run differently. To get a universal Yukawa coupling, higgsinos and gauginos should be 
included. However, we will be satisfied with the present approximation in this paper 
for two reasons. First, the Yukawa coupling ht is determined by the top mass which is 
not known very precisely at the moment. Second and even more important, the most 
significant effects of ht appear in conjunction with the soft squark mass parameters 
mQ^mfj^m^^ (see below). Since these are unknown, there is a large uncertainty in the 
calculation in any case. Once the squark masses have been measured and the top mass 
is known more precisely, the gauge coupling ht{ji) should be fixed at 1-loop level in 
terms of the top pole mass. 

We will work in the MS scheme (with the scale parameter Ji) and take Tr 1 = 4. For 
the squark and quark loops included in vacuum renormalization the results agree with 
those in the DR-scheme, often used in supersymmetric theories. 



4 Dimensional reduction 

Let us first recall the expansion parameters of dimensional reduction |TH|. Since all 



the calculations are IR-safe, no non- analytic powers of masses can appear. In fact, the 
expansion proceeds just in powers of 



as at zero temperature. We include only the quarks and squarks of the third generation 
in the loops affecting vacuum renormalization, whereas the corrections e.g. from gauge 
bosons, suppressed by (^^/(IGTr^), are neglected. 

In addition to expanding in coupling constants, we make a high-temperature expan- 
sion in the mass parameters. This requires that the soft supersymmetry breaking mass 
parameters satisfy 

ml, ml, mj^, ml, ml,, < (27rTc)^ (4.2) 

The limit (|4.2| ) implies that m^<2'7rTc, so that the results of this paper cannot be 
directly continued to the limit mA oo studied in ^, ^. Actually, the limitation 
on mA is not as important as that on the squark mass parameters, since the latter are 
associated with larger coupling constants. 

To keep track of the validity of the high-temperature expansion, we will at some 
points display also the leading correction terms. The critical temperature is Tc ~ 100 
GeV so that we shall assume mg, m(/, m^, ?^a^300 GeV. We also assume that the 
masses generated at the electroweak phase transition as well as the masses associated 
with possible colour and charge breaking minima are small compared with 27tTc. 
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If some of the soft masses are large, one cannot use the high-temperature expansion. 
Instead, one should evaluate the corresponding integrals numerically. If m ^ SvrTc, one 
can also integrate out these degrees of freedom in the sense of the vacuum decoupling 
theorem [p7| . 

The actual dimensional reduction proceeds by writing down the general form of the 
effective 3d theory and then determining the 3d coupling constants by matching the 
Green's functions in the original theory and in the 3d theory. The degrees of freedom 
of the effective theory are the bosonic degrees of freedom of the original theory. The 
temporal components of gauge fields become Higgs fields in the adjoint representation. 
The structure of the 3d theory is determined by gauge invariance. The 1-loop calcula- 



tions needed are a straightforward application of the rules in |T^. We just write down 
the graphs and the results below. 

Since the complete bosonic sector of MSSM is rather large, we display only the part 
interacting with the SU(2) and Higgs degrees of freedom explicitly. We recall that after 
trivial rescaling with T, the dimension of bosonic fields in 3d is GeV^/^ and that of 
the couplings gl,g'i,h'^3,hl^ is GeV. At some points, we denote new parameters with 
the same symbols as the old ones, to avoid increasingly cumbersome notation. Higher- 
order operators suppressed by the temperature and coupling constants are neglected. 
Finally, let us recall some basic notation: 

cb = ln(47r) - 7s ^ 1.953808, = cs - 2 In 2 ^ 0.567514, (4.3) 
= 21n^-2cB, Ly(/i) = 21n|-2cs. (4.4) 

The 3d effective theory consists of the following parts: 

1. The temporal components of the original gauge fields become scalar fields in the 
adjoint representation, and the spatial components remain gauge fields. Relevant for 
the present discussion is the part 

W = ^Ft^Ft^ + ^{D,A^oY + Irn^AoA^o + ^(^.C'o^)' + Irnl^foC^, (4-5) 

where D^Aq = diAQ + g^e'^^'^A'^AQ. When only quarks and squarks are included, the 3d 
fields are related to the renormalized 4d fields by 



1 



9' fr , 1 



IGtt^ V ' 2 



(4.6) 
(4.7) 



where mass corrections suppressed by tuq/ {2'kTc)'^ were neglected. 

The gauge coupling can be most easily obtained from the graphs (qqqq), (SS), (SSSl), 
(SSS2), (SSSS) in Fig. |l|.c. Here also a redefinition of the Higgs fields, given in (|4.12|) - 
( [4.13| ), is needed. After the redefinition one gets 



gl = g\n)T 



9' (r , 1 



167r2V ' 2 



(4i 
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(P) 



(PP) 



X-/ / '-^ r 

(b) >c > . :x 

(PP) (PPP) (PPPP) 



(PP) (PPPl) (PPP2) (PPPP) 



Figure 1: The generic types of graphs needed for dimensional reduction of (a) wave 
functions and masses, (b) scalar couplings and (c) the gauge coupling. Wiggly lines 
are vector propagators and dashed lines represent generic propagators of particle type 
P=Q, U, D, S, H, A, C, g, q, f. Here Q, U, D denote the corresponding squarks, 
S is a squark in general, H is a Higgs doublet, A and C are the SU(2) and SU(3) 
gauge fields, g is a ghost, q a third generation quark and f a general fermion. For the 
coupling constants, 1-loop dimensional reduction is directly related to 1-loop vacuum 
renormalization and hence only squarks and quarks are considered in the internal lines. 
For the masses, the thermal screening terms proportional to are not related to 
vacuum renormalization and hence we include all the modes with m<27rT in the loops. 



new) 



so that gl{A\A])^ 

The values of '^Ao5'"^c'o known 0], but for completeness we write them 

here as well. These terms contain only screening parts not related to vacuum renor- 
malization, so that we include the complete spectrum of the model in the loops. With 
the notation in Fig. |l|.a, the graphs contributing to are (ff), (gg), (AA), (HH), 
(QQ), (A), (H), (Q); to mco contribute (ff), (gg), (CC), (SS), (C), (S). For illustration, 
the leading mass terms are also shown: 



m 



Co 



9l 



2 Up = 

3 3 



1 + 



Up 



+ 



— ]T^ + 

6 y 



87r2 
1 

8^' 



[m-^ + 1712 + 3"^q) 



+ ml + 2ml) 



(4.9) 
(4.10) 



In eqs. ( [4.9|) , ( [4.10|) , n^? = 3 is the number of fermion families and Ng is the number of 
scalar doublets interacting with the gauge fields in question. 
2. The quadratic terms of the Higgs sector are 



^Higgs 



{D,Hy{D,H') + {D,Hy{D,H^) + i\{D,H')\D,H^) + H.c 



+ mlH'^H' + mlH^^H^ + 



mUH'^H' + H'^H'] 



(4.11) 
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where Di = di — ig3T°-A'^/2 and is the 3d gauge couphng. The graphs contributing to 
the 2-point Higgs correlators are (qq), (AH), (SS), (A), (H), (S) in Fig. [I|.a. According 
to the general strategy of including only quark and squark loops in terms related to 
vacuum renormalization, we neglect terms multiplied by g'^ everywhere except for the 
screening parts proportional to T^. The new fields are then related to the renormalized 
4d fields in the MS scheme by 

C(3) dl + wl- 



1 + 



T 



167r2 
1 



3hiLfm) + 



l + -—3hiLfm) + 



1287r4 T2 
C(3) u 



s ^ 



167r 



1287r4 T2 



The mass parameters, on the other hand, are 



(4.12) 
(4.13) 



o(now) 

mf ' 



9 (new) 



m 



2 (new) 



12 



+ 



+ 



+ 



+ 



+ 



3 



167r2 



-hi + -q' 



rp2 



C(3) 



1287r4T2 



3dl{ml 



(^: 



(4.14) 



167r2 



3.2 1 ' 
-hi + -q' 
4* 4^ 



C(3) 



1287r4T2 



o 2/ 2 



+ m^) 



• 2 I 2\ 2 



(4.15) 



m 



12 



167r2 



-{h^ + hl)ml2Lf{fi) + {usWs - dses)Lb{fi) 



C(3) 



1287r4T2 



— 3dses(m% + 



D) 



1 



(m^ + ^2 + + 



12 



(4.16) 



where we have shown terms up to quadratic order in the masses. From the high- 
temperature expansion, one would also get terms of the form K^rn^ in addition to the 
'u277i2-terms shown above, multiplying the coefficient C(3)/(1287r^T2), but these terms 
have been neglected. In fact, the renormalization structure of the theory suggests the 
parametric convention m ~ /i^T, Ug ~ hiT , according to which all the terms involving 
C(3)/(1287r^) would be of higher order. Nevertheless, we keep the terms shown since 
numerically the mixing parameters might be larger than some of the masses. 

The Lhifi), Lf{fi)-teTms on the first rows of eqs. (|4.14|) - (|4.16|) cancel the running of 
rriKfi), m2(/i), ml2{p') so that the 3d masses are RG-invariant at 1-loop order. More 
precisely, the effect of the Lb{fl), Lf{p,)-teTms is to run the MS mass parameters to a 
certain scale (It, which need not be the same for all the parameters. For instance, if 
there are only bosonic contributions, then it can be seen from eq. (4^) that /It ~ 7T. In 
Sec. ^ the running parameters rnKfl), mKfl), 'ml2{p') are expressed in terms of physical 
parameters and fl so that the /i-dependence cancels in the 3d parameters. 
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Finally, the parameter in eq. (|4.11|) is 

C(3) 



1287r4T2 



(4.17) 



If the (yf^-corrections from the Higgs fields were included, there would also be a term 
proportional to g'^m\2 inside the parentheses in (f4.17|) . 
3. The quadratic terms needed in the squark sector are 



(4.18) 



The graphs contributing to the 2-point correlator (Q^'Q) are (AQ), (CQ), (SH), (A), 
(C), (S), (H) in Fig. |l].a; for {U*U) and {D*D) the interactions with SU(2) gauge fields 
are missing. The fields in ( ^.18| ) are related to the original fields by 



( D* _D)(°*^™) 



1 

f 



1 . 2r / N C(3) ul + dl + wl + el 



167r2 
1 

IGtt^ 
1 



167r2 



4:glLb{fi) + 



3847r4 T2 
C(3) ul + w! 
1927r4 T2 
C(3) rf^ + e^- 



1927r4 T2 



(4.19) 
(4.20) 
(4.21) 



where the g'^-teims have been neglected. The terms proportional to g^ represent the 
contributions within the present theory, and are due to gluon loops. 
The mass parameters in eq. ( [4.18|) are 



m 



U3 



mm 



+ 



^Q(^^) - T^^bifi) ht{ml + ml) + hl{ml + m^) 



167r2 

+ < + dl + el + wl-^ 



2 2 



C(3) 



1287r4T2 



2/ 2, 2\ij2/ 2, 2\i 2/ 2, 2\| 2/ 2, 2\ 

M^(m2 + mj;) + d,(mi + m^) + w;,(mi + m^) + e,(m2 + m^) 
1 



+ 2{usWs - dses)m\2 - -{ut + di + ei + wi)m^ 



2\_2 

Q 



Y^LtifL) [2hl{ml + ml) + 2(^2 + ^2) _ _^2^2^ 



C(3) 



1287r4T2 



2ul{ml + ml)+2wl{ml + ml) 



1 



"^d(^) - Ye^^'''^^) [2/ib(mi + mg) + 2(4 + ej - -^^m^ 
2 , C(3) 



1287r4T2 



2d2(mf + m^) + 2e2(m2 + m^) 



4rf,e,mi2- -(ci^ + e^)m^ 



(4.22) 



(4.23) 



(4.24) 
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The Lfe(/i)-terms cancel the running of the soft masses mQ{fi),mfj{p,),m'jj{p,) so that 
the 3d masses, hke the Higgs sector masses, are RG-invariant at 1-loop level. It should 
be noted that the running may be noticeable; for instance, the parameter mfj which 
may be small has a running proportional to hfrriQ, so that the relative effect may be 
significant. As in the case of Higgs mass parameters, one should hence renormalize the 
squark mass sector at the 1-loop level to remove the //-dependence (1-loop corrections 



to the stop mass have been calculated in |28|). However, since the squark masses at zero 



temperature are not known at the moment, we will not perform any renormalization in 
the present investigation. Instead, the parameters mg^fiT), ^ui.P'T), ^'oifiT) produced 
by the Lft(/2)-terms in eqs. (|4.22| )- (|4.24| ) are replaced by the tree-level values. 



4. The interactions of Aq and Cq with the Higgs fields are produced by the graphs 
(qqqq), (SS), (SSSl), (SSS2), (SSSS) in Fig. |l].c. There are terms of the form 

L*^*^^ = hiA'^A'^H^^H^ + h2AlAlH^^H^ (4.25) 

existing already at the tree-level, and terms generated radiatively, 

= h^AlAl{H^^ + H^^ H^) (4.26) 

The coefficients related to Aq are 

^ 4 1 287r4 T ' ^ ' 

3 ^\(3) UsWs-dsCs , . 

- "4128^ T • ^^-^^^ 

Note that in the Standard Model extra terms of the type g'^h^T / (IGvr^) are generated in 
hi, h2 through quark loops |T^, but in MSSM these terms are cancelled by the squark 
loops. If the (/^-corrections from Higgs fields were included, there would be a term 
proportional to g'^m\2 in eq. ( |4.29| ). 

As to the coefficients Ci,C2, the quark loops (qqqq) give the contributions 

5ci = -^A', ^<^2 = ~^,9lhl (4.30) 

as in the Standard Model, but these are cancelled by the squark loops (SS), (SSSl) in 
Fig. |I].c, as for hi, h2. Hence there only remain the small terms 

^1 = -^^T^' (4.31) 

- ~ 64vr4 T ' (^-^^^ 
= _^|C(3)^.^.-4e 

^ 647r4 T ^ ' 
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5. For the quartic self- interactions of the Higgs fields, the most general gauge- 
invariant two Higgs doublet potential |^ is generated at the dimensional reduction 
step. Since we assumed all the parameters to be real in the original Lagrangian, the 
potential is somewhat simplified, being of the form^ 

V = \i{H^^H^f + \2{H^^H^f + X^H'^H^H^^H^ + X^H^^H^H^^H^ (4.34) 
+ X^{H^^H^H^^H^ + H.c.) + XqH^^H\H^^H^ + H.c.) + XrH^^ H\H^^ H'' + H.c). 

To give the expressions for Ai, . . . , A7, we use the functions 

f[ 1 1 

•'Pb 



f2{ma,mb) 



(p2 + m^)(p2 + m^) {p 



2\2 



(4.35) 



C(3) ml + ml 



2^2 



+ 



m 



167r2(27rT)4 



f^{ma,mb,m^) 



% (p2 + ml){p'^ + ml){p'^ + ml 



(4.36) 



C(3) 



C(5) ml + ml + ml 



f4{ma,mb,mc,md) 



10247r6 



+ 



m 



167r2(27rr)f 



^6 (p2 + ml){p'^ + m2)(p2 _^ ml){p'^ + ml) 



(4.37) 



C(5) 



10247r6T4 



+ 



m 



167r2(27rT)f 



where the sum-integral is over the non-zero bosonic Matsubara frequencies in the MS- 
scheme. In the numerical computations we keep only the constant part in the function 
/a, to be consistent with the fact that terms of the same parametric form come from 
the redefinition of fields and the higher mass contributions were there neglected. The 
part subtracted in the definition of /2(?7iaj^fe) is 



1 



Pb 



,2^2 



167r2 



1 



+ Lbifi) 



(4.38) 



The graphs needed are (qqqq), (SS), (SSS), (SSSS) in Fig. |T|.b. In addition, the 
redefinitions of fields according to eqs. ([4.12|) -( fi.l3|) give contributions. After the 
redefinition, the parameters are (^| = gl/T): 



Ai ^ ^ i 1 

T 8 8 167r2 



Lbm)-Lfm)]{--g'-3ht + -g'hl) 



(4.39) 



Y^9^f2{mQ,mQ) - ^hf,[f2{mQ,mQ) + f2{mD,mD)] + ^g^ hi f2{mQ , mq) 



^We recall that the identity {H^^ H^){H^^ H^) + [H^^ H'^){H^^ H^) = {H^^ H^){H^'< H^) reduces 
the number of independent combinations. 
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3 3 1 

■g'^dlf3{mQ,mQ,mD)- -dif^{mQ,mQ,mD,mD) - -w^f4{mQ,mQ,mu,mu)>, 



7 8 8 1 67r2 



L,{-^)-Ls{-^)]{--g'-^ht + -g'ht) 



(4.40) 



.2(^,2 

3 



C(3) 



4^ ^"- + ^^^1287r4r2 



^ i ~T^^^-^2(mQ, mcj) - -/i^ [/2(mQ, mg) + /2(m[/, mu)] + jg'^hlf2{mQ, mq) 
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fsimu, mu, mq) + h{mQ, mq, mu)\ - -g^lhimq, mq, mu) 



As 
T 



3 3 3 1 

+ ^g^(^1f3{mq,mq,mD)- -ujf4{mq,mq,mu,mu) - -e^Ji{mq,mq,mD,mD)>, 

(4.41) 



4 4 167r2 



L,(/.) - L^(/.)] [--g' - 6hX + -g\hl + hi) 



-^g^f2{mq,mq) 



- Sh'^fhl [f2{mq, mq) + f2(mu, mo)] + |/(^t + hl)f2(mq, mq) 

+ 3htwlf3{mu, mu, mq) + Sh^dlfsimq, mq, mo) + Shle^Jsimo, mo, mq) 

+ 3hlulf3{mq, mq, mu) + Qhthiugds + esWs)h{mq, mu, mo) 

+ ^^^(^s - ul)h{mQ, mq, mu) + ^/(e^ - dDhi^q, mq, mo) 

- 3{usds + esWs)'^f4{mq, mq, mu, mo) - 2>ulwlf4{mq, mq, mu, mu) 

- 3d1e1fi{mq,mq,mD,mD)y 



gl 



T 2 ^ 167r2 



L,{-^) - L;(/x)] ^-g' + Qhlhl - \g\hl + hi 



(4.42) 



+ \g\< + + + 1281^^2 + \^^9*h{mq, mq) 

+ ^hlhl [f2{mq, mq) + f2{mu, mo)] - ^d^iht + hl)f2{mq, mq) 

+ Shfwlfsimq, mq, mu) - Sh^d^fsimq, mq, mo) + Shle^Js{mq, mq, mo) 
- 2>hlulf^{mq, mq, mu) - Qhthiugds + esWs)h{mq, mu, mo) 

+ ^9\ul - w^^)f3{mQ, mq, mu) + ^g^i^l - e^Mmq, mq, mo) 

\2 



+ 3{usds + egWsf fiimq, mq, mu, mo) - ^ulwlf^imq, mq, mu, mu) 
- ^dlelf4{mq,mq,mD,mD)Y 

Y ^ [^-\ulwlfi{mq,mq,mu,mu) -^dlelfi{mq,mq,mD,mD)Y 



(4.43) 
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T 



fsimD, iriD, itiq) + fz{mQ, uiq, uid) 



(4.44) 



T 



(4.45) 



3 3 

+ ^g'^UsWsfsiruQ, ruQ, ruu) + ^g'^dsCshimQ, tuq, tud) 

- 3Uswlfi{mQ, niQ, mu, rriu) + 3esdlfi{mQ, mg, mo, mc] 

3h^UsWs[f3{mu,mu,mQ) + f3{mQ,mQ,mu] 

3 3 

- ^g'^UsWsfsi'mQ, rriQ, mu) - -g'^dsCsfsimQ, uiq, uid) 

- SwsulUiniQ, niQ, mu, mu) + 3dselU{mQ, mg, m^, m^,) 



The terms in the curly brackets will be useful in Sec. |^ as well, which is why they have 
been separated. 

6. Cubic interactions of Higgs fields and squarks are of the same form as in the 
original theory: 

i^cubic = UssH^^QaUa + d^sH^^QaD^ + CssH^^QaD^ + WssH^^QaU^ + H.c. . (4.46) 

Since the terms Us,ds,Ws,es are unknown, it is not so important at the moment to 
calculate the 1-loop corrections to the tree-level formulas. Just as an illustration of the 
structure that appears, let us give the 1-loop terms proportional to ht, within the 
present theory: 



Ms3 



ds3 



es3 



Ws3 



Usifi) - 


1 

167r2 




3 

+ hthbds)Lb{fi) + -h'fusLf{fi) 


5 




(4.47) 


dsifi) - 


1 

167r2 


{Qhlds 


3 1 
+ hthbUs)Lh{ft) + -hldsLf{fi) 






(4.48) 


esifi) - 


1 r 

167r2 i 


{Shies 


3 

- hlcs + hthbWs)Lb{fl) + -hie 


sLf{ji) 




(4.49) 


Wsifi) - 


1 

167r2 


{Shfw^ 


3 

, - hlws + hthbes)Lb{ji) + -hi 


WsLf{n) . 


(4.50) 



T 

7. Quartic interactions of Higgs fields and squarks are at tree-level of the form 

-Z^quartic = -^glH"^* H'J'Q*^^Qia{25il8jk — Sij6kl) 

+ hUH'^H^KU^ + H^^Qo^QiH') + hUH'^H'DlD^ + H^^Q^QIh') 
+ hMH^^H^UlD^ + l^.c). (4.51) 

In principle it would be important to calculate the 1-loop corrections especially to ht^ 
since it affects the transition quite significantly. However, as stated above, this is not 
accessible within the present framework, since vacuum renormalization of the top quark 
mass cannot be used to simultaneously fix the /it's appearing in different places in the 
Lagrangian beyond tree-level. Moreover, the top mass is not known very accurately. 
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and the effect of hts comes togetfier witli mu^ whicfi is not known at all. Hence we take 
the couplings here only at tree- level. Then the couplings are as written in eq. ([4.51 ) 
with 

hl = h',T, hl, = hlT. (4.52) 

8. Finally, there are many terms not interacting directly with the SU(2) gauge fields 
and Higgs fields. We will not show them explicitly, since they enter the further inte- 
grations only at 2-loop level. Nevertheless, in some cases the higher-order corrections 
are important; a particularly relevant example [0 is discussed in Sec. 0. We just fix 
one more notation here: the strong coupling constant in 3d is denoted by g^g^ and is 
gl, = gUfiT)T. 



5 Integrating out squarks, Aq and Cq 

The bosonic theory discussed in Sec. ^is still rather complicated, although simpler than 
the original theory. However, generically many of the fields appearing are massive at 
the phase transition point. Such fields can be integrated out in 3d. It appears that 
in some part of the parameter space, all the squarks together with the adjoint scalar 
fields Aq, Cq can be integrated out. 

More specifically, the requirements for the integration to be valid are the following. 
First, the phase transition should be weak enough so that the neglected higher-order 
operators are not important. Second, the perturbative expansion for the parameters 
of the effective theory should converge. The first requirement should be reasonably 
well satisfied when a;>0.03 which is the region we are studying. Let us investigate the 
second requirement in some more detail. 

Integrating out Aq, Cq gives roughly the expansion parameters 

9l 9h 



From eqs. ( [4. 91) , ( [4.10|) one sees that these are small numbers, below 0.05 (to be more 
precise, the expansion parameter of C,^-integration might be slightly larger due to 
colour factors, but on the other hand g"^^ appears first only at 2-loop level). With the 
trilinear couplings (which have the dimension GeV^^^ in 3d) are associated expansion 
parameters of the type 



Ws3 



(5.2) 



which are very small for small mixing. The largest and most important expansion pa- 
rameters are related to the strongly interacting squarks. There the expansion proceeds 



in powers of (see Sec. |T0|) 



Js^^J!k_ (5.3) 

and correspondingly for the other squarks. Roughly, the factor 4 in the denominator 
of (|5.1| ) is compensated in (|5.3| ) by colour factors. The terms in eq. ( ^.3|) are of order 
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0.3 if 



^U3 



(5.4) 



in which case the neglected 2-loop terms are expected to give a correction of about 
20% (1-loop corrections may sometimes be almost as large as tree-level terms). In the 
present Section we shall assume that mc/>50 GeV so that the expansion in ( ^.3|) should 
still be useful (for the other squarks, we assume mg ~ ~ 300 GeV). The case of 



smaller mu is discussed in Sec. 10 



In the case that all the squarks and the A^- and C^-fields can be integrated out, the 
new theory will be 



L = -Ft^Ft^ + {D,Hy{D,H') + {D,Hy{D,H')+^[{D,Hy{D,H') + }l.c.] (5.5) 

+ X^iH'^H^H^^H^ + H.c.) + X^H^^H\H^^H^ + H.c.) + XrH^^H^H'^H^ + H.c). 

Although the notation for the parameters is the same as before, the parameters have 
changed from the previous theory. 

The graphs needed for calculating the parameters have a simple relation to the graphs 
needed in the dimensional reduction step. The quark contributions do not exist any 
more. The squark graphs remain precisely the same. In addition, there are the extra 
graphs with Aq, Cq in the internal lines, of the same type as for squarks but without 
cubic interactions with the Higgs fields. 

From the graphs (SS), (AqAq) in Fig. [^.a, one gets that the new fields are related to 
the previous ones by 



/ j,N (new) 



1 + 



(new) 



{a^a 



9l 



+ 



9l 



167rmQ3 247rm^o 



2 

HI 



'H2j 



where 



^Hl 



^H2 



1 + 



1 + 



i. 



W 



s3 



+ 



r/2 

"s3 



All V (mQ3 + mf/s)^ (mgs + itld^ 



.(. 



u 



s3 



+ 



The parameter ^ is changed to be 

1 / Us3Ws3 43653 



(5.6) 

(5.7) 
(5.8) 

(5.9) 
(5.10) 



^(new) 



Att V (mQ3 + m(73)3 (mQ3 + tuds) 



(5.11) 
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The new gauge coupling can be derived from the graphs {AqAq), {AoAqAqI), (SS), 
(SSSl), (SSS2), (SSSS) in Fig. |T[c and is 



2(new) 
^3 



9l 



9l 



9l 



167rmQ3 247rm^o 

SO that 5^3^ ^{A^AA = g^iA^AA. The new mass parameters are given by 



^2(new)^2 



HI 



mi 



4tt 



Air 



w 



s3 



+ 



m. 



2(now) ^2 



'H2 



m 
3 



2 - -^^^niAo 



Air 



An\mQ3 + mu3 mQ3 + uids 



htsiniQs + mu3) 



u 



s3 



m 



2(new) r7 ry 

12 ^H\^H2 



4tt 



+ 



^s3 



Us3Ws3 



ds3Gs3 



(5.12) 



(5.13) 



(5.14) 



(5.15) 



Arc V mgg + mf/g TTlQs + mD3 ' 

For the scalar coupling constants, one can to a large extent use the results in 
eqs. ( |4.39| )-( p:.45| ). The squark graphs and the combinatorial factors are precisely the 
same, but the integration measure and the parameters appearing have changed. The 
fermion graphs are missing, but the Aq, Co-graphs have to be included. Hence the 
graphs are (SS), (SSS), (SSSS), (AqAq), (CqCq) in Fig. g.b. We will display only the 
part arising from Aq, explicitly; the rest can be read from eqs. ([4.39|) -( ^.45D and is 
indicated by the curly brackets below. The replacements to be made in (|4.39| )-( p:.45| ) 
are that the functions f2, f3, fi are replaced with those defined in eqs. (|5.16| )- (|5.18|) 
below; g, ht, Us, . . . ^ 93, hts, Uss, ■ ■ •; and m Q,my , mp- ^ mgg, nius, mns- 

The integrals appearing, analogously to ( [4.35|) -( p:.37|) , are 

dp- 



f3{ma,mi„mc) 

1 



dp 



(p2 + m^)(p2 + m^) 47T{ma + rrib) 
1 

(p2 -I- m1){p'^ + ml){p'^ + m^) 
1 



ATT{ma + mb) {ma + mc){mb + rric)' 

^ ijp' + m^) {p^ + m^) ijP' + m1)(jP' + 
ma + mb + mc + rrid 



(5.16) 



(5.17) 



Att {ma + mfc) {ma + (m^ + m^) {m^ + mc) (m^ + m^) (m^ + m^) 
The integration measure here is 

d'^p 



dp 



{2Txy 



d = 3-2e. 



(5.18) 



(5.19) 



18 



There is no divergence in f2{fna,'mb) in 3d, so that it was not necessary to subtract 
anything in the definition in contrary to the 4d case. 
With the notation introduced, the new parameters are 





= Ai 


- ^hlf2imAo,mAo) - 


8clf2{mco,mco) 


+ {/2- 


-/4}, 


(5.20) 


X (new) r74 
A2 ^fj2 


= A2 


- 3hlf2{mAo,mAo) - 


8clf2imco,mco) 


+ {/2- 




(5.21) 


^ (new) y2 y2 


= A3 


- Qhih2f2{mAo,mAo) 


- 16ciC2/2(mco, 


^Co) + 


{/2 - h 


}, (5.22) 


^ (new) ry2 ry2 


= A4 


- 6/13/2 (m^o'"^Ao) - 


16c3/2(mco,mco 


) + {/2 




(5.23) 


■V (new) r72 rv2 

^H1^H2 


= As 


- 3hlf2{mAo,mAa) - 


Sclf2{mco,mco) 


+ {/2- 




(5.24) 


, (new) r73 ry 

'^e ^Hl^H2 


= Ae 


- 6/11/13/2(^^0,^^0) 


- 16ciC3/2(mco, 


^Co) + 


{/2 - U 


}, (5.25) 


, (new) ry ryS 

A-j ZjhI^H2 


= A7 


- 6h2hf2{mAo,mAo) 


- 16C2C3/2(mCo, 


^Co) + 


{/2 - U 


}. (5.26) 



Here we have displayed the terms arising from field redefinitions on the LHS of the 
formulas. The factors Zhi, Zh2 are given in eqs. (|5.9|), ( p. 101 ). 



6 Diagonalization of the two Higgs doublet model 



The theory in eq. (|5.5|) can still be simplified. The phase transition should take place 
close to the point where the mass matrix has a zero eigenvalue. Then generically the 
other mass is heavy. Recall that at tree-level the sum of the eigenvalues of the mass 
matrix is ml+m^ = m^, and at finite temperature one gets positive thermal corrections 
to the masses. Hence one may integrate out the heavier Higgs doublet as well. In order 
to do so, we first diagonalize the two Higgs doublet model. 

We make the diagonalization in two steps. In the first part we rotate and rescale the 
fields so that the term 

e[(Ai^')^(Ai^') + H.c.] (6.1) 

disappears from the Lagrangian in eq. ( ^.5|) . In the second part we rotate the resulting 
fields so that the non-diagonal mass term 



m 



12 



{H^^ + H^^ H') (6.2) 



disappears. Then the resulting theory will be 



L = -i^^.i^^. + (A0)^(A0) + (A^)^(A^) + m|0V + ^; 



1 

4" 

+ Ai(0V)^ + A2(^^^)^ + A30V^^^ + K(t)^ee^(t) 

+ A5(0^^0"f0 + H.c.) + A60V(0"^^ + H.c.) + A70"f0(0^0 + H.c.). (6.3) 



It should be noted that for small values of the squark mixing parameters, ^ in ( |6.1| ) is 
very small so that the first part of the diagonalization is numerically inessential. 
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The first part of diagonalization proceeds by writing 





= J_ '(^l _ ^)-l/2^1(new) _)_ (^1 _|_ ^^-l/2^2(new)" ^ 


(6.4) 




= J_ — ^)-l/2_f/'l(now) (1 ^)-l/2_f^2(now)" _ 


(6.5) 


Expressed in terms of the new fields, the term in eq. (|6.1|) vanishes. The other param- 


eters become 






2(new) 

nil 


= (l-0"'K + m2-2my/2, 


(6.6) 


2(new) 


= {l + i)-\ml + ml + 2ml^)/2, 


(6.7) 


2(new) 
"^12 


= {l-er"\ml-ml)/2^ 


(6.8) 


■I (now) 
^1 


= (1 - 0"'(Ai + A2 + A3 + A4 + 2A5 - 2\ - 2A7)/4, 


(6.9) 


(new) 
^2 


= (1 + 0"'(Ai + A2 + A3 + A4 + 2A5 + 2A6 + 2A7)/4, 


(6.10) 


(now) 
^3 


= (l-e')"'(Ai + A2 + A3-A4-2A5)/2, 


(6.11) 


■I (now) 
A4 


= (l-a~'(Ai + A2-A3 + A4-2A5)/2, 


(6.12) 


-1 (now) 
^5 


= (l-a-i(Ai + A2-A3-A4 + 2A5)/4, 


(6.13) 


(now) 
^6 


= (1 - 0-^(1 - e)-"\)^i - A2 - Ae + A7)/2, 


(6.14) 


■1 (now) 
A7 


= (1 + 0-^(1 - eV\>^i - A2 + Ae - Ar)/2. 


(6.15) 



In the second part of the diagonahzation, we write 



cosa + sina 9, 
— sina (h + cosa 9. 



The angle a is chosen so that 
tan2a; = 



2ml^ 



rrio — m 



2 ' 
1 



sin2a 



2ml^ 



ml 



m^)^ + 4m| 



(6.16) 
(6.17) 



(6.18) 
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It should be reiterated that at 1-loop level in the MS scheme the 3d mass parameters 
are finite, so that we need not worry about renormalization at this point. As a result 
of the rotation in eqs. (|6.16| ), ( |6.17|) , the action is of the form in ( |6.3|) . The new mass 
parameters, obtained from those in (|6.6|)-([H^), are 



mr, 



rrii + 1712 ~ V ("^1 + 4m' 



12 



ml + ml + 



ml 



(6.19) 
(6.20) 



Abbreviating cos^a 



c^, sin^a 



s^,cos2q; = c2,sin2Q; = s2, the matrix M giving 
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the couplings as A'-"'^^-* 



MX from those in ([OD -dOSD, is 



M 



s2V4 
s2V4 



c s 

s2V2 s2V2 

s2V2 s2V2 

s2V4 s2V4 -s2V4 

s2c^ -s2s^ -s2c2/2 

s2s2 -s2c2 s2c2/2 



c* + s^ 
-s2V2 



s2V4 

s2V4 

-s2V2 

c^ + s^ 

-s2V4 

-s2c2/2 

s2c2/2 



32^2 
32^2 
-s2^ 


-ch2 


-sh2 


sH2 


cH2 


s2c2 


-s2c2 


-s2^ 


s2c2 


-s2c2 


c^ + s^ 


s2c2/2 
c2 - s22 


-s2c2/2 
s22 - 


-s2c2 


s2c2 


s2^ - 


C2 - S22 



(6.21; 



7 Integrating out the heavy Higgs doublet 



In eq. ( 6.18|) the angle a has been chosen such that the field is light at the phase 
transition point, as can be seen from ( |6.19| ). Then the heavy field 9 can be integrated 
out. The expansion parameter is 

(7.1) 



9l 
ATTrrif, 



which is very small in the cases we are studying (recall that mg'>,mAiT). It should be 
noted that arises for the first time at 2-loop level, whereas at 1-loop level only the 
scalar self-couplings appear. 

When 9 is removed, the resulting theory is just the 3d SU(2)+IIiggs theory: 



1 



^3d 



F^F^ + (A</>)^(A</') + m20t0 + A^(0t^ 



(7.2) 



For this theory there are non-perturbative lattice results available, so that one need 
not go any further with perturbative methods. 

Since the interactions in the starting point, eq. (|6.3| ), involve vertices of the type 
and 6'0^, there are non-standard graphs needed in the construction of the effective 
theory. Numerically these graphs may not be very important since the relevant coupling 
constants Ag, A7 are not large and are suppressed by the large mass me in the results. 
Nevertheless, conceptually the way to include Ag, A7 has to be addressed. It should be 
noted that while in eqs. (|4.44|) - (|4.45|) Ae, A7 are much smaller than Ai, . . . , A4 for small 
mixing parameters, in general this is no longer true in the theory of eq. ( |6.3| ) due to 
the redefinitions of fields in Sec. 

Let us start with the wave function normalizations. The wave function (j) does not get 
normalized in the integration, since at 1-loop level there are no momentum-dependent 
contributions to the 2-point correlator (00) from the heavy modes 9. Due to the graph 
in Fig. |I|.a, the wave function A" becomes 



b\ (new) 



The gauge coupling is changed to 



[a^a] 



1 + 



9l 



ASnrng 



2(new) 
^3 



-V 



(7.3) 



(7.4) 



21 



(a) 



(b) ^xcix; 



{0) 



(d) 



Figure 2: The graphs needed for integrating out the heavy Higgs doublet from the 
3d two Higgs doublet model. The solid line represents the heavy field 9. Graph (a) is 
a contribution to the mass parameter m^, graphs (b) are contributions to the scalar 
self-couphng A^, (c) is an induced 6-point function, and (d) is a mixing term generated 
at 1-loop level. 

For the scalar mass parameter, the diagram in Fig. |^.a gives 

m5(-) = mJ-^(2A3 + A4). (7.5) 

The scalar coupling constant receives contributions from the diagrams in Fig. ^b to 
become 

A<^ = Ai - (A^ + A3A4 + \\l + 2\l + I2\l - I2\h) . (7.6) 

The couphng A2 would enter only at 2-loop level. Let us discuss the result for the 
coupling constant in some more detail. 

The contributions involving A3, A4, A5 come from graphs of the type [QQ^ in Fig. ^.b, 
and are standard. 

The contribution proportional to Ag comes from the graph (^0), involving a light 
field in the internal line. In principle, one might think that Ag contributes at order Ag 
only to the 6-point function depicted in Fig. ^.c. Such a contribution, however, has a 
momentum-dependence: 

C6 ~ (7.7) 
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To construct a local effective theory, one wants to expand in the momenta. Naively 
one might think that it is justified to expand eq. (|7. 7| ) everywhere in p^/m^ since the 
effective theory only involves the mass scale -C mg. This naive procedure is wrong, 
since if the ^-propagator is expanded before integration in the graph {6(f)) of Fig. ^.b, 
one only gets the suppressed contribution 

Xl'^. (7.8) 

In reality, the dominant contribution of (6'0) is of order Xl/rrig. Only when this larger 
contribution is explicitly included in the reduction step by the graph {6(j)), can one 
expand in the momenta in eq. ( [7.7| ). Then, in fact, the 6-point operator 0^ can be 
neglected in the effective theory, since it only leads to contributions suppressed by 
m^/mg < 1. 

The phenomenon explained is of course the same which appears in the dimensional 



reduction step at 2-loop level [0, |T^] when one is comparing a naive integration over 
non-zero Matsubara frequencies with a matching procedure for the construction of an 
effective theory. The former method, containing exclusively heavy modes in the internal 
lines, leads to a non-local theory. The latter method leads to a local effective theory, 
but light fields have to be included in the internal lines of some graphs. In the present 
case, the difference appears already at 1-loop level. In principle, a systematic way to 
account for these effects is to split the light fields into low- momentum modes 0|p|<A 
and high-momentum modes 0|p|>a; then only the heavy fields and the high-momentum 
modes of the light fields need be included in the internal lines. 

The contribution (6) in Fig. |^.b is even more exotic than (6'0). It cannot even be 
generated from an effective potential for the 0-field alone as the other contributions, 
since the graph is reducible. This contribution arises because the vertex involving A7 
induces a mixing between (p and 6 at 1-loop level, as shown in Fig. ^.d. This mixing 
does not vanish in the limit that mg is large, but grows as X^mg. Nevertheless, it is still 
possible to construct order by order an effective theory of the type in (fT^), containing 
the light fields only and giving the same light Green's functions as the original theory. 

In the configuration {9) in Fig. ^b, the induced mixing contributes to the 4-point 
function of the 0-fields at the same order of magnitude as {00), (Ocj)). To reproduce 
this contribution in the theory of eq. (|7^ ), the graph (0) has to be included in the 
reduction step. Working at 1-loop order, one may expand the momentum dependence 
of this graph, but going to 2-loop order, the graph obtained from (0) by contracting 
the rightmost light field with one of the other light fields has to be included in the 
calculation of the mass parameter to order A^. 

Finally, let us recall from Sec. that the parameter relevant for baryogenesis in the 
theory of eq. (|7.2| ) is the dimensionless ratio x = A^/^ff at the phase transition point. 
The temperature dependence of X^/T,gl/T is weak: at the dimensional reduction step 
the dependence comes only through logarithmic 1-loop corrections. In the heavy scale 



integrations a larger dependence is induced since e.g. h'^^/mu^ ~ h'^T / ^ mfj + jT^ 
depends on T. We estimate the critical temperature from the condition 

ml = 0, (7.9) 
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which gives sufficient accuracy for the present purpose. In particular, note that im- 
posing eq. ( |7.9| ) in the 3d effective theory generally involves the next-to-leading correc- 
tions |30| to Tc in terms of the original coupling constants, arising from the heavy 3d 
modes. The numerical results obtained for x are discussed in Sec. ^. 



8 Vacuum renormalization 

To complete the program of dimensional reduction, one has to fix the running param- 
eters appearing in Sec. ^ in terms of zero temperature pole masses and cross sections. 
We reiterate ffist the general strategy adopted in the present paper. 

1. The Higgs mass parameters (/i), m^/i), '^i2(/u) are determined in this Section. 

2. The running of mg(/i), m^(/i), m|i(/i) can be read from eqs. ( [4.22D -( p:.24D . Since 



the squark masses are not known at present, we do not fix the running parameters in 
terms of pole masses here. However, once the squark masses have been measured, the 
renormalization must be properly performed, since especially the possibly small stop 
mass parameter roy has a significant effect on the phase transition. In the present 
work we use tree-level values for mQ(/iT), m^(/ir), m|,(/iT). 

3. For the present type of investigations, the running gauge coupling is most conve- 
niently fixed in terms of the muon lifetime |T^. In the complete MSSM, the result could 



be extracted from a calculation of the type in [^. However, as stated before, within the 
present theory the gauge coupling appearing is not universal beyond tree-level if other 
than squark and quark loops are included in vacuum renormalization. Hence there is 
no use in going into elaborate investigations; we will rather fix 5f(/2 = 200GeV) = 2/3 
and include only the running due to quarks and squarks. 

4. The U(l) gauge coupling is taken at tree-level, and is fixed to be g' = 1/3. 
At 1-loop level no difference is made between g"^ and g"^ = g"^ + g'^- Since vacuum 
renormalization is related to the neutral sector of the theory, we will use the numerical 
value of g^ in the loops calculated with the SU(2) interactions. 

5. The Yukawa couplings ht{fi), hb{fi) can in principle be fixed in terms of the top 
and bottom pole masses. However, within the approximations of the present paper, 
there is no universal hflfi) at 1-loop level. Hence we fix also ht,hh at tree-level. For 
fixed tan/3 = ^2/^1 and fi, we take 

q TUf 1 , q rrih 1 

ht = ^ — —^. hb=^ — 8.1 

V 2 mz sm p ^2 mz cos p 

6. The mixing parameters are also running parameters. However, they are not very 
important for the phase transition, at least if small, and they are not known. We fix 
them too at the tree-level through 

Ws = -fiht, (8.2) 

Cs = fJ^hb, (8.3) 

Us = -htAt-Ws cot (3, (8.4) 

ds = —hbAjj + estanp. (8.5) 



24 



It general, Wg and are arbitrary soft supersymmetry breaking parameters and hence 
fi may be interpreted as something different from the supersymmetric mass parameter 
in the superpotentiaL With the conventions in 
in the broken phase, 



Ml 



m 
m 



ui 

2 

U12 



m 
m 



1/12 
2 

U2 



]^), the squark mass matrices 

(8.6) 



and analogously for A4jj, are given by 



m 



ui 



m 



tL 



rri' 



2 I 2,-^2 na 2 



m 



m 



U12 



1 / 



— TJltA 



m 



Dl 



1 



+ ml 



m 



D12 



--^(esf2 - dsVi] 



-m| cos2/? . 
rribAb. 



m 



D2 



tR 



rrij 



mjj + m; 



1 1 



(8.7) 
(8.8) 
(8.9) 
(8.10) 



We next concentrate on fixing mKfi), m2(/2), 'n^^2(/^) terms of the pole masses m^, 
uiA, Tnz of the lightest CP-even Higgs particle h, the CP-odd Higgs particle Ha and 
the Z-boson Z^, respectively. Going to the classical broken minimum determined by 
rriKfi), niKp,), m^gl/^) the scalar couplings, one can calculate the tree-level masses 
'^lifi), m^ifi), m|(/x): 



Here 



-Tn\{fi) + [mKfi) - m^(/i)] /cos2/3 (/i) 
1 



f8.111 



m^(/i) +m^(/i) 



sin2/3 (/i) 



2m?2(/i) 



4m^(/i)m|(/i) cos22/5 (yu) 



f8.12) 



mf(/i) -hmi(/i)' 

Adding to the tree-level expressions the 1-loop self-energies Ilh{k'^,ft) = {h{k)h{—k)) , 
IlA{k'^,fi) = {hA{k)hA{—k)) , and the transverse part Ilz{k'^,jj') of the Z-boson self- 
energy {Z^{k)Z^{—k)) , evaluated at the corresponding poles, gives the physical fi- 
independent masses: 



rrir 



m\{ji) - Y[A{-m\,Ji), 



(8.13) 



From these equations one can solve for ml{fj,),m\{fi),m%{fj,) for given m\,m 
Using the expressions inverse to ( ^.11 ), 

^lil^) = - +"^l(/i)]cos2/3(/i)}. 
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^2(/^) = ^{mi(/i) + [mi(/i) + m|(/i)]cos2/3 (/.)}, (8.14) 



2 

^i2(/^) — — m^(/i) sin2/5 (/i) 



where 



sin2/?(^) = — pT\/[mi(/i) - ml{n)][ml{fi) - mUfi)], (8.15) 

one then gets the desired expressions for m\{ji),m\{ji),m\2{ji). Note that at tree- level 
the running parameters equal the physical parameters and eq. (|8.15| ) implies the known 



relation < min(myi,m2), but at 1-loop level the running parameters are no longer 
pole masses and consequently the Higgs pole mass can be considerably larger than 
implied by the tree-level bound. 



A few comments are in order. First, in solving eqs. ( |8.13| ) it is technically convenient 
to keep tan/? fixed rather than the Higgs mass. The parameter tan/3 is not a physical 
quantity, though, and depends on the gauge, on the scheme and on Ji. When tan/3 is 
fixed, TTih is not an input parameter any more but comes out as a result from (|8.13|) 
for given tan/3 and Ji. 

Second, we choose to use the physical pole masses as the mass parameters in the 
self-energies n(/c^,/i). This reduces the higher loop /i-dependence of the result. In this 
procedure, one also needs the unknown pole mass of the heavier CP-even Higgs mass 
niH (in the tadpole diagrams). The corresponding self-energy can be trivially obtained 
from YihihP', p,), see Appendix A. Then one has to add the additional unknown uih and 
the additional equation 

= ^Hifi) -'^H{-mH,JJ') (8.16) 

to the three equations ( ^.13| ). This system of four equations [fixed: /i, tan/3, niA, mz; 
unknown: m/j, mn, m\{ji), m^^Jj)] eqs: (|8.13| ), (|8.16| )] is easily solved by iteration. 

Third, for the present purpose it is sufficient to work strictly at 1-loop level. Numer- 
ically, it is important in some regions of the parameter space to include higher-order 
corrections, arising for instance from the mixing induced at 1-loop level between the 



tree-level mass eigenstates h and H ||3^ 



The 1-loop self-energies have been calculated in the literature [0, |2^ in detail, but 
for completeness we also display the formulas used here in Appendix A. Using the tree- 
level expressions in terms of m^, m\2, mq, mfj, m|, for the parameters appearing, 
it is straightforward to verify explicitly that the //-dependences produced for mf, m^. 



mfg through (|8.13|) , (|8.14| ) agree to leading order with the ones in ([4.14|) -( p:.16|) . It 



should be noted, however, that numerically the remaining higher order /i-dependence 
may be as important as the leading order one. Fortunately, the /i-dependence of the 
coupling constants determining the Higgs mass bound is smaller than that of the mass 
parameters, see below. 



26 



9 Numerical results 



Combining the results for vacuum renormalization with the formulas for dimensional 
reduction and heavy scale integrations one can study the values of x in the phenomeno- 
logically allowed part of the MSSM parameter space. The phenomenological constraints 
on the squark sector, relevant for the present analysis, have been discussed in 0, ||, ^ . 
First, there exist lower bounds on the masses of the weakly interacting squarks. This 
gives a lower bound on itlq, the stronger constraint arising from sbottom. More im- 
portant, the relative mass splitting 

(m?^-m?J/m?^ (9.1) 

of the left-handed squarks is constrained by the parameter Ap(t, h) not to be too 
large. Since m?^ contains and m? contains and both contain mg, this also 
acts as a lower limit for mg. We will take mq = 300 GeV which should satisfy the 
phenomenological constraints for our reference value = 175 GeV within the accuracy 
of the present calculation. At the same time the chosen value of ttlq still lies within 
the applicability of the high-temperature expansion. 

For the right-handed stop mass parameter mfj there appear to be no phenomeno- 
logical lower bounds apart from the absence of charge and colour breaking ||^. On the 
other hand, one cannot take too small values within the applicability of the integra- 
tions in Sec. |^, since then the expansion parameters in eq. ( p.3| ) grow large. We take 
mu = 100 GeV as a reference value. The mixing parameters are taken to be zero at 
the reference point. At = Ab = fi = (we stress again that fi in (^.2|)-( p3D is not really 
the supersymmetric mass parameter affecting the chargino and neutralino masses, and 
hence a small value for it is acceptable). 

We fix the renormalization scale used in vacuum renormalization and dimensional 
reduction to be /2 = 200 GeV. The dependence of physical quantities (such as x) on 
fi is formally of higher order than the accuracy of the present calculation. In practice 
there is some dependence which may be used to estimate the accuracy of the results, 
see below. This dependence arises for instance since hf^ depends on fi, having been 
fixed at tree-level. 

We will next vary the CP-odd Higgs mass between 50 and 300 GeV and in- 
spect the values obtained for x, for different values of the lightest CP-even pole Higgs 
mass rrih- In particular, the dependence on mt,mu and the mixing parameters At, 
around the reference point is of interest. 

In Fig. ^the value of x is shown as a function of for the reference set of parameters 
with three values of tan/3 (thin lines) and three values of the lightest Higgs mass 
rrih (thick lines). We recall from eq. (|2.2|) that the requirement for a strong enough 
phase transition to sufficiently suppress the sphaleron rate in the broken phase is x < 
0.03 — 0.04. First, we notice that the best region for baryogenesis is a heavy CP-odd 
Higgs particle, as is already known Second, for the reference parameters, even the 
region x < 0.03 can be reached with a sufficiently small Higgs mass (although then 
our results are less reliable, see Sec. ^and below). This should be contrasted with the 
situation in the Standard Model where it appears that no Higgs mass is possible [jl9 
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Figure 3: The effect of the CP-even and CP-odd Higgs masses and mA on x [x is 
defined in eq. (|2.2| )). All the numbers are in GeV. The mixing parameters have been 
set to zero. With the thin lines, we show as an alternative parametrization the value 
of function of tan/? (at = 200 GeV) in the present scheme. 



Hence the situation has definitely improved in the MSSM. However, the Higgs mass 
needed in the MSSM would be rather small, m/j<70 GeV. This might soon be excluded 
experimentally. 

In Fig. ^ the effect of the top mass is shown, for fixed tan/3 (thin lines) and fixed 
rrt/i = 70 GeV (thick lines). In contrast to the Standard Model (see Fig. 27 in |I9| ) 
a large top mass makes the situation more favourable for baryogenesis. The reason 
for the difference is that in the MSSM the top Yukawa coupling also appears in the 
dimensionally reduced 3d effective theory through squark interactions. The corrections 
induced for the scalar self-coupling are large and negative as seen in eqs. ( [4.40| )-( p:.42| ). 
However, since the corrections are large, they are also sensitive to the precise value of 
hfi- A way to estimate the reliability of the results is their /i-dependence, which exists 
since hfi is fixed only at tree-level. By varying Ji from 200 Gev to 300 GeV for fixed 
mhi the change in x is less than 3% for mt < 175 GeV. For rrit = 190 GeV the change 
is about 10%. Hence the calculation becomes less reliable for large top mass. 

It should be noted that we have kept ttlq fixed when varying rrit. In fact, if rrit is 
larger, then also mq is likely to be larger, in order to keep the mass difference of left- 
handed stops and sbottoms small as required by phenomenological constraints [|, 
This effect would compensate for the increase in the strength of the transition with 
rrit 0- III principle, rrit also directly affects the running of mQ{fi),m'lj{p,) to /i = fir, 
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Figure 4: The effect of the top mass on x. Here the top mass is taken at tree-level, 
in accordance with the other uncertainties in the calculation. The thick lines are for 
constant mh = 70 GeV, the thin lines for constant tan/5 = 2.0. The mixing parameters 
have been set to zero. 

but these effects have been neglected here. 

In Fig. 1^ the effect of varying the squark mass parameter my is shown, again sep- 
arately for fixed tan/3 and for a fixed Higgs mass. A smaller mjj makes the situation 
more favourable, as was already noted in [§, ^, ^j. In even negative values for m\j 
were considered. Here we cannot go to that region since then the squarks are not heavy 
any more and the effective theory is different, see Sec. 0. Nevertheless, one can see 
how the effect starts to arise. It is seen that for irtu = 50 GeV even a Higgs mass in 
the region rrih ~ 75 GeV seems possible. 

Finally, in Fig. ^ the effect of the mixing parameters is presented. Comparing with 
Fig. ^, one can see that /i has very little effect {x is just slightly reduced at tua ~ 300 
GeV). Indeed, for large rriA the mixing is determined exclusively by the combination 
At appearing in the squark mass matrix. The effect of At is that a large value makes x 
larger. Phenomenologically, this is somewhat unfortunate ^ since one might wish to 
have a non-zero mixing in order to get smaller squark masses, which might help with 
the i?b-problem. The sbottom mixing parameter Ab has practically no effect at all. 

The overall conclusion is that for mh<mw and small mixing, the transition might be 
strong enough for ■mu<50 GeV and mA>200 GeVQ. If the mixing is larger, one would 

^ In it was proposed that another favourable region is at small , independent of tan/3. From 
Fig. it can be seen that x is indeed almost independent of tan/3 for ~ 50 GeV (this feature 
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Figure 5: The effect of mu on x for constant rrth (thick hnes) and tan/3 (thin lines). 
Note that at tree-level, the right-handed stop mass m^^ at zero temperature is given 
by m? 



nij 



+ ml for vanishing At and g' , see ( ^^ -( ^78|) . 



need an even smaller mjj. These results agree to large extent with [4-7]. In any case, 
the situation has certainly improved with respect to the Standard ModelQ. 



10 The effective 3d theory in the case of light right- 
handed stops 

It has been stressed in 0, that small values of are phenomenologically allowed 
and are favourable for electroweak baryogenesis. On the other hand, it was pointed 
out in Sec. ^ that, if m\j is small, one cannot integrate out the right-handed squarks 
and the relevant effective 3d theory is not the simple SU(2)-|-Higgs theory discussed 
in the previous Sections. In this Section we discuss in more detail at which point the 
3d integration is no longer reliable and what the relevant 3d theory is then. We also 



persists also for values of tan/3 larger than shown in Fig. g). However, x is still much too large and 
rn/i much too small, as can be seen from the m/j — 60 GeV curve. Hence the effect proposed in 
does not take place close to our reference point. 

^ It should also be noted than in the Standard Model there is a critical Higgs mass above which 



the phase transition ceases to be of first order |33 . In the MSSM, on the contrary, there exists an 



upper bound on to/j, and in some cases (e.g. in the vicinity of our reference point) all possible Higgs 
masses result in a first order transition. 
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Figure 6: The effect of the mixing parameters on x for constant mh (thick hnes) and 
tan/? (thin hnes). The notation At in the figure stands for Af. If /i = or if ttla is 
large so that fi has httle effect, the results are symmetric under At —A^ so that only 
positive values are shown. For smaller ttia, the increase in x is smallest when the signs 
of and At are the same. 

indicate some non-perturbative effects which might arise beyond the perturbative ones 
discussed in P, 0. 

In 0, the case mu^, ~ was investigated. From eq. ( |5.3| ) it is clear that then 
the integration does not work at all. To get a quantitative estimate of the m^/s still 
allowed, one should compare 1-loop and 2-loop contributions to different parameters 
of the effective theory: /i^a-corrections first arise at 1-loop level and (y's'a-corrections at 
2-loop level, so that a comparison of tree-level and 1-loop results does not reveal much. 
The dominant 2-loop contributions were identified in 0. 

The dominant 2-loop effect [0 is due to graphs of the type 

V21 = -2glVssv{mi^,mi^,0) (10.1) 

in the notation of 0, |23|. Here essentially m?^ = m^g + h'^v'^/^ (in [0 the coupling 
constant appearing in this formula is h^sin"^ P due to the hmit 00). According 

to eq. (82) of ||15|, the term in (|10.1|) affects the dimensional reduction step of Sec. |^ 
only by changing the mass parameter by terms of the type 

^^l-^29lh'tln^. (10.2) 
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These terms are related to the running of hf{ft) in the 1-loop thermal correction and 
hence their inclusion requires 1-loop renormalization of the top quark mass. In any 
case, these terms only affect the critical temperature and thus are not very important. 
The (finite) 3d-part of the 2-loop contribution, on the other hand, is |^ 



29lVlUm,^,m,^,0) = S^^fm?^ ^ f In + (10.3) 



Now, if mu3 is large enough and the transition is not too strong, this term can be 
expanded in powers of /i^f|/(2m^3). The first term, proportional to V2, changes the 
mass parameter m| at 2-loop level, and is not very important. The second term is 
quartic in f 2 and changes the coupling A2 by 

For clarity, we have here kept the coupling constants in their 4d normalizations so that 
powers of T are written explicitly. The change in eq. ( |10.4| ) is negative, reducing the 
coupling constant A2 and consequently making baryogenesis more likely. Hence, as 
long as the expansion converges, this 2-loop correction works in a favourable direction 
also in the framework of the effective SU(2)+Higgs theory discussed in Sees. §-0. 
However, when the effect becomes stronger, the convergence becomes worse and the 
higher order operators generated become important. The right-handed stops can no 
longer be integrated out but act as light degrees of freedom. 

The expansion parameter of [/-field integration can be estimated by comparing the 
2-loop term in eq. (|10.4|) with the corresponding 1-loop term in eq. ( ^.21|) : 



o /, 47-12 
16n mu3 



Hence the expansion parameter is roughly 

6Xf 2 glT 



[10.6) 



Consequently, to get convergence one needs mu^^T. 

What would be the effective theory if mu^<T and the integration does not converge? 
Let us assume that itla —>■ 00, the squark mixing parameters are small and rriQ is 
relatively large as required by phenomenological constraints for a realistic top mass. 
Then all the other squark degrees of freedom apart from U can be integrated out in 
the dimensionally reduced 3d theory. What remains can be written down immediately 
using 3d gauge invariance: 

T ^ rpa rpa , ^ y^yAyryA 

+ {D^H)\D^H) + mlH^H + XuiH^Hf 
+ {Dl U)\Dt U) + ml U^U + \u {U^Uf 
+ 'j^H^HU^U. (10.7) 
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Here Df = di — ig3T^A1/2 and Df = di — igsz^^Cf /2 (in this effective theory we 
have denoted the complex conjugate of the original [/-field by U). At tree-level 73 = 
hi sm^(3T, Xu = c/|T/6 and ml = ml + (%|/9 + h^/Q + hj sin^ 13 / Q)T'^ . 

The steps needed for a more precise derivation of the theory in eq. ( |10.7| ) are in 
principle the following. First, make dimensional reduction as in Sec. ^but in the theory 
where -^00. In a precise study, it would be important to consistently include 
all the 1-loop corrections to the Yukawa couplings ht appearing in different places. 
Second, integrate out the heavy fields Q, -D, Aq, Cq as in Sec. H. Finally, make vacuum 
renormalization in order to fix the MS-parameters in terms of physical parameters. In 
particular, one should renormalize and ht^Jj) in addition to the Higgs sector 

parameters by calculating the stop and top masses at 1-loop level. All these steps are 
straightforward and parallel the ones presented in Sees. ^j-p. 

Let us stress that perturbatively the theory in eq. ( |10.7| ) reproduces the 1- and 2- 
loop results making the dominant effects in P, 0]. In fact, the 3d theory also contains 
a resummation of IR-safe higher-loop contributions, so that it is expected to be more 
precise than direct perturbative calculations in 4d. More important, eq. ( |10.7| ) contains 
all the IR-problems of the theory and could be used for 3d Monte Carlo simulations. 
No such simulations are available at the moment for the complete theory. However, 
one can try to use the knowledge obtained from simulations of the SU(2)+Higgs sector 
to get some insight into the properties of the complete theory. We make two guesses. 

1. In the simulations of the 3d SU (2) -(-Higgs theory it was found that in the sym- 
metric phase the relevant degrees of freedom are non-perturbative bound states ||T^, |2^, 



22| , |3J]. The mass of e.g. the scalar bound state may differ much from the perturbative 
value, let alone from the tree-level value. If ml is positive at Tc so that the SU(3)-part 
of eq. ( p.0.7| ) is in its symmetric phase, one might expect the same phenomenon to take 
place here, only the effects would be stronger than in the SU(2)-sector. This is impor- 
tant since the results of |^ strongly depend on ml and assume a tree-level value 
for it. In particular, the non-perturbative mass might be significantly larger than the 
perturbative and tree-level masses, in which case the light degrees of freedom of the 
theory at the phase transition point might again be described by the 3d SU(2)+Higgs 
model. This time, however, the derivation of the effective theory would have to be 
non-perturbative. 

2. The symmetric structure of eq. ( [10. 7| ) opens other interesting possibilities. At 
the phase transition, m\j is close to zero, and if is also rather close to zero as 
proposed in 0, one might end up in a situation where also the charged and coloured 
field U acquires a non- vanishing expectation value at some point during the transition. 
This kind of a multi-stage transition might naturally alter the mechanism of baryoge- 
nesis. Requiring the absence of colour and charge breaking during the transition, some 
constraints on the parameters were given in ^ . 

A precise investigation of the possibilities proposed will have to wait for a detailed 
perturbative derivation and a lattice investigation of the theory in eq. ( p.0.7| ), as well 
as for experimental data on the values of the unknown parameters. 
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11 Conclusions 



We have constructed super-renormalizable 3d effective field theories describing the 
thermodynamics of the electroweak phase transition in MSSM. The derivation of these 
theories is perturbative and free of IR problems. The effective theories can then be 
used for further perturbative investigations if IR problems are beheved under control, 
or better still, for non-perturbative Monte Carlo studies. 

It was found that in a part of the parameter space, it is possible to reduce the 
effective theory to a 3d SU(2)+IIiggs theory for which there already exist lattice results. 
However, it generically appears that when the reduction can be done that far, the 
transition tends to get rather weak for realistic Higgs masses. Pushing the parameters 
into the region of a stronger transition (a smaller right-handed stop mass parameter 
mfj), the convergence of the 3d heavy scale integrations gets worse. 

It hence seems that for a strongly first-order transition, the relevant effective 3d the- 
ory may be more complicated than SU(2)+Higgs. A particularly appealing possibility 
is a model containing an SU(2) scalar doublet and an SU(3) scalar triplet. If mfj indeed 
turns out to be small, this effective theory should probably be studied in more detail. 
As far as the derivation of the theory is concerned, the most important pieces missing 
at the moment are the expressions for the parameters /if(/iT) and m^(/iy) in terms of 
zero-temperature physical parameters beyond tree-level. The calculations required are 
straightforward and parallel the calculations presented in the present paper. 
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Appendix A 

In this appendix we give the formulas used for vacuum renormalization in Sec. H More 



complete expressions can be found e.g. in [p2| , |28| (for compact approximation schemes 



and some 2- loop corrections, see e.g. and references therein). 

The calculation of the 1-loop self-energies is organized as follows. We first shift 
the fields to the classical broken minimum. Then the mass eigenstates are identified 
and the 1-loop graphs needed for {h{k)h{—k)), {hA{k)hA{—k)) and {Z^{k)Z^{—k)) are 
calculated. In particular, the tadpole graphs have to be included since we are not at 
the exact quantum minimum. According to the general strategy of this paper, only 
quarks and squarks of the third generation are included in the loops. 

For fixed tan/3, the location of the classical broken minimum of eq. (|3.3| ) is obtained 
from 

2mz cos P ^niz sin /3 / . n 

vi = z , V2 = z . (A.l) 

9 9 
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At the broken minimum, the mass eigenstates corresponding to the physical neutral 
Higgs fields H, are obtained from the fields in eq. (|3.2|) with the rotations 



hi 



cosa h + sina H, 
— sina h + cosa H, 
cos/5 — sin/3 Ha, 
sin/? hco + cos/? Ha- 
At tree-level the angle a here is given by 



hi 
hi 
hi 



sin2a 



sin2/? 



2 2"' 

- mi 



cos2a = cos2/? 



2 2 
— 

ml -ml' 



(A.2) 
(A.3) 
(A.4) 
(A.5) 



(A.6) 



At 1-loop level we use the physical pole masses for m^, m|, ml, m\ and the angle a 
is determined from the expression for cos2q; . After the redefinitions (|A.2|) - (|A.5|) , the 
graphs contributing to {hh), {hAhjC) and [Z^Zy) can easily be identified. 

The formulas arising are slightly complicated by the fact that the left- and right- 
handed squarks mix. For the mass eigenstates, we will use the notation m^)^ 
defined by 

(A.7) 



m 



u± 



m 



Ul ^U2 

are in (|8.7|) - 
1 



and correspondingly for where m\ji, 

1 



2 + 4m^ 



r/12 



.10). We also denote 



mfj^ - ml_ 



= — —■ (A.8) 

Some standard integrals often appearing are denoted as follows. For \mi — m2\ < 

k < mi + m2, 



FH{k;mi,m2) 



m? — , m-i m? + , mi 
In — H k i In — 



fc2 



m2 



2 2 
mf — m2 



m2 



{mi + m2y — k'^Jk'^ — {mi — m2y arctan 



(A.9) 

\Jk'^ — {mi — m2Y 



{mi + m2) 



Especially, 



Fff(mi;m2,m2) = 2 



2^4^ 



1 arctan 



(A.IO) 



where r = m2/mi and r > 1/2. Outside the displayed kinematic region, an analytic 
continuation is needed, and from that we only use the real part in calculating the 
masses. The imaginary parts arising are small. We also define a function Fz arising in 
the calculation of the Z-boson self-energy: 

^ ^ r\ 2/ 2 2 2 \ 1 TTXl r\2/2 2 2\i 

2mi[mi — m2 + m^) in h 2m2[m2 — mi + m^) m 



Fz{mi,m2) 



V2m\ 



mz 



+ 



m 



I - 2m\{ml + ml) + {m\ - mlf\ \l + In 



m% 



mz 

m? -|- mi mi 
m — 



+ FH{mz;mi,m2) 



H — mr 
3 ' 



mim2 

3m|(mi + ml) — {ml — m2)^ 



mf — ml m2 



(A.ll) 
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(a) 



(qq) 



(SS) 



V / 



(S) 



(b) 



(qq) 



(SSI) 



/ 9 



(SI) 



(S2) 



(SS2) 



h?V2 ^ h?V2 

(SS3) 



(SS4) 



htV2 



Figure 7: The graphs contributing to (a) the pole mass mz of the Z-boson and (b) 
the masses mh,mH,mA of the CP-even Higgs bosons and the CP-odd Higgs boson. 
In the internal lines, solid lines are quarks and dashed lines are squarks of the third 
generation. Examples of couplings appearing are also shown. 



The irreducible graphs needed are shown in Fig. and the tadpole graphs in Fig. ^. 
We leave out the common factor l/(167r^) in the formulas below. 

The contributions to {h{k)h{—k))\j^2^_^2 from the graphs in Fig. |^.b are 



(qq) : 



(SSI) : 



(SS2) 



— sin a 



(6m? - mV) In — ^ + 2m? + (4m? - ml)FH(mh] rrit, rrit) 
mf 



+ — > hb, "^fe, sin^a — > cos^a |, 

3(ut sin^a + w"^ cos^a — UsWg sin2a ) < In ■ 



2 \2 
U2) 



1 + 6u{mlj^ + mfj^) In ^^^^ + ^^(m/,; mu-, mu+) 



+ {u 
3 



d., w. — > e,, sin a ^ 



cos^a , sin2a — sin2a ,U D^, 



16 



g (^2 sin a + cos a + V1V2 sin2a ) In 



ml 



+ —g sin a + Vi cos a + V1V2 sin2a )6u< (mr/2 — mrj 
16 



In- 



mr. 



mt 



+FH{mh; mu-, mu-) + (m^2 - ^u+)'^ 



In 



mr. 



(A.12) 



(A.13) 



m 



FH{mh;mu+,mu+] 



+2m' 



U12 



1 + ln 



mr 



u+ 
mu- 



+{u^d}, 



■ + Suimfj^ + mfj^) In + Fnimh] mu-,mu+) 
mu-mu+ mu+ 



(A.14) 
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(SS3) : Qh^V2 sin a <\n 



mu-mu+ 



(SS4) : 



I { U T, 2 • 2 2 2 

+ 1 /if — > /ife, f 2 sm a — > ti]^ cos a 
— ^g^h^ sin^a + -f if 2 sm2Q; ) In 



U^d}, (A.15) 



— -g'^h^[ylsm'^a + -fif2 sin2a; 



2 ^2 , ™4 

+ "'(712 



In 



rrir. 



m 



u- 



+FH{mh;mu-,mu- 



+ 



[m 



U2 



In 



u+ 



+FH{mh; mu+, 'rnu+) \ + hh, sin^a cos^a - 



U^d}, (A. 16) 



(SI): 



cos2a mfr, (in + 1 
4 ^^V m| 



3 

— -^^ cos2a 6u 
-{U^D}, 



[m 



U2 



mlj_)m'lj_ In 



mz 



mz 



(A.17) 



(S2) : ^h; sin^a {mjj^ + mjj^) In + 1 



nir 



— o/ij sm a \rnu_ m h m 

+ \ht — ^ ^fe, sin^a — > cosset , — >• D j. 



(A.18) 



The contributions to {hA{k)hA{—k))\ 



from the graphs in Fig. |^.b are 



(qq) : -^Kj cos^/? {2mi - m\) In + 2mi - m^Fnim^; m) 



mi 



+ \ht — > /ifo, nit nih, cos^(3 — i> sm^(3 |, 



(A.19) 



(SSI) : 3(u^ cos^/5 + wl sin^P — UgWs sin2/3 ) 



In- 



mu-mu+ 



+ 1 + 6u{m1j^ + mfj^) In h FuimA, ruu-, mu+) 

mu+ 



(SI) 



+ 1^5 ^ dsi Ws Cs, cos^/9 ^ sin^/3 , sin2/5 
3 



sm2p,U ^ d], (A.20) 



cos2(3 m'fj^ ( In + 1 



mi 



3 ~ r 
+ -g'^ cos2/? (5c/ [(m^2 - ml-)T^u- In 

(S2) : ^hl cos^/? (m^i + m^2) fin + l) 



mz 



(m^2 - mjj+)mjj+ In 



777^ 

(A.21) 
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—Qhf cos^B (mfr In — - — h mfr. In 

The contributions to the transverse part of {Z^{k)Z^{—k))\ 
Fig. |^.a are 



A;2 = — m| 



(A.22) 

from the graphs in 



(qq) 



(SS) 



r,2 



(3m^ - m|) In + + (m^ - mz)FH{mz; rrit, rrit] 



mt 3 



1 _2 / 2 ^ 2 \ 1 



(A.23) 



ml 



+2m^i2^z(^c/-, mu+)] + {U d}, (A.24) 



2 ^2 ,^2 ^2,21^!^ 



[m^2 - "^;7_)"^c/_ In 



m 



C/2 "''U+)"''U+- 



mz 



(A.25) 



To the contributions in eqs. ( |A.12| )-( |A.25| ) one has to add the tadpole contribu- 
tions from Fig. |], since we are working around the classical minimum. The tadpole 
contributions to n/i(/c^, /i), nyi(A;^, /i), n2(A;^, /i) are 



n 



(tad) 



^ COs2a ^(tad) 



COS2/5 



Vi sin^a + V2 sin2a /2 



m 



H 



viCh + ( 



V2 cos^a + f 1 sin2a /2 



m 



H 



n 



(tad) 



f 1 cos^a + f 2 sin2a /2 f 1 sin^a — f 2 sin2a; /2 



--^2cos2/3 

f2 sin^a + f 1 sin2a /2 ^ f2 cos^a — vi sin2a; /2 



ml 



V2Sh 



(A.26) 

V2Sh , 

(A.27) 



n 



(tad) 



Vi cos^a — V2 sin2a /2 ?;i sin^a + V2 sin2a /2 



viCh 



+ 



t>2 sin^a — t>i sin2a /2 112 cos^a + vi sin2a /2 



m 



V2Sh 



(A.28) 



Here 



2hfm^ — — + {dsCs — UsWs) cot/5 



h]{mui' + mc72^) + ^g^{mui' - m^i') 



In^V + l 

mi 
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/ \ / N / \ 

\ J \ / \ / 

\ht I I hlv2 I g^V2 

Figure 8: The tadpole graphs needed in vacuum renormahzation. The closed loops 
contain quarks and squarks of the third generation, and the single line may contain 
either of the CP-even Higgs particles. Representative coupling constants are shown. 

10A2 2i . „,2/ 2 1 , 2 1 ^U- 

— lln^m^m \rbn^\mjj,m VmTj_m 

niz ^ vfiz rrtz 

'mu+ 2 1 'f^u- 



+ GSuUsiug + Ws cot/3 ) (m^T, In — — — mij_ In 

^ mz mz 

+ Qdo^si^s — dsCotl3){m%,,hi — ^ — m?5_ In — ^ 

^ mz mz 

3=2 



3 



2^ 



2 2 \ 2 1 ^U- / 2 2 \ 2 1 

U2 - m^J)m^_ In {m^^ - m^^)m^^ In 

mz mz 

mo- ,2 2x2, 



+ 7;fSD {ml, - ml_)ml^ In ^ - (m^ - ml^)ml^ In ^ , (A.29) 



, „_ D+ "''D+ 

mz mz 



Ch = \^ht —>■ hb, mt — > mb, Ug ^ dg, Wg ^ — e^, cot/3 — * tan/3 , f/ ^ d|. (A. 30) 

The contributions to the heavier CP-even Higgs mass are obtained from the contri- 
butions to the lighter CP-even Higgs mass by changing m| m|^, sin2a; — sin2a , 
cos2a — cos2a , sin^a cos^a inside the 1-loop formulas for Ilh{—m1, fi), as can 
be seen from eqs. ( |A.2| ), ( [A.3| ). 
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Erratum (February, 1999) 

In Appendix A, some loop contributions were erroneously omitted from the correlator 
{h{k)h{—k))\,2^_^2, Eqs. (A.12)-(A.18). These contributions, which do not have any 

h 

scale dependence, are proportional to squark mixing parameters and are thus very small 
for the small mixings considered. Nevertheless, they should in principle be included. 
In the notation of Fig. 7.b, the loops omitted are of the form 

Us ^ ^ q'^vo Us ^ ^ h'jv2 

(SS5) (SS6) 
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The results for these contributions, to be inserted between Eqs.(A.16),(A.17), are: 



3 r 
(SS5) : — y=g^{—UsSma + WsCOsa){viCOsa + V2sma)5urnlj^2\^^ 



(SS6) 



+5u 



U2 



m 



ui 



+(m^2 - niu_)FH(mh; mu-,mu-) + {m^u2 - TT^u+)FH{mh, mu+, mu+) 
+(m^i - rrvfj2)FH{mh; mu-, mu+] 



+ |— tis sina + Wg cosa — >• dg cosa + sina, — >• d|, 
12htmt{—UsSm^a -\- -WsSm2a)5um^-^2 



(A.ie.a) 



X 



f21n^^^ + FH{mh;mu-,mu-) - FH{mh;mu+,mu+) 

hf,, rut — > ruh, Us sin^a — > —ds cos^a ,Ws ^ es,U — > (A.16.b) 

Numerically, these graphs only affect Fig. 6 where non- vanishing mixing was consid- 
ered. The curves for fj, — 200 GeV do not change, since the omitted contributions are 
proportional to At, Ai,. The largest effect is for At = 100 GeV: even then, the thin line 
(fixed tan/?) remains essentially unchanged. The thick line (fixed m^), comes down by 
~ 0.005. All the conclusions remain unchanged. 

I thank M. Losada for bringing these omissions into my attention. 
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